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The  use  of  solar  energy  for  heating  greenhouses  in  the  cold  period  for  a  Mediterranean  climate  is  an 
important  issue.  In  this  paper,  the  thermal  performance  of  a  new  solar  air  heater  using  a  packed  bed  of 
spherical  capsules  with  the  latent  heat  storage  system  in  east-west  oriented  greenhouse,  is  analyzed  and 
discussed.  The  excess  heat  in  the  greenhouse  was  stored  in  the  packed  bed  through  the  diurnal  period  and 
extracted  at  night.  An  experimental  comparative  study  was  conducted  in  two  greenhouses  installed  in  the 
Research  and  Technologies  Centre  of  Energy  (CRTEn)  in  Tunisia.  The  greenhouse  heat  balance  in  different 
components  of  the  greenhouse  with  a  heating  system  (cover,  canopy,  soil,  inside  air  and  packed  bed  solar 
collector),  has  been  used  to  investigate  the  impact  of  the  phase  change  material  (PCM)  on  the  greenhouse 
temperature  and  humidity.  Results  show  that  the  nocturnal  temperature  inside  the  greenhouse  equipped 
by  a  heating  system  exceeds  the  temperature  inside  the  conventional  greenhouse  by  5  °C.  The  relative 
humidity  was  found  to  be  of  an  average  10-20%  lower  at  night  time  inside  the  heated  greenhouse. 
The  nighttime  recovered  heat  of  the  solar  system  attained  31%  of  the  total  requirements  of  heating.  An 
economic  evaluation  was  made  to  find  the  payback  period  of  the  new  solar  air  heater,  the  potential  in 
saving  fuel  energy  and  the  reducing  carbon  dioxide  emissions. 
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1.  Introduction 

Greenhouses  have  been  used  in  various  forms  of  protection  for 
the  plants  when  it  would  have  been  impossible  to  grow  outside 
because  of  the  adverse  climate,  diseases  and  pests.  The  greenhouse 
permits  to  achieve  maximum  returns  from  cultivation  by  creating 
better  growing  conditions.  The  growth  factors  namely  light,  humid¬ 
ity,  temperature  and  nutrition  should  be  delivered  and  maintained  at 
optimal  levels  [1,2].  To  keep  the  optimum  growing  environment, 
greenhouses  should  allow  high  light  transmittance,  sufficient  venti¬ 
lation  efficiency,  low  construction  costs  and  low  heat  consumption 
[1,3  .  With  the  present  increase  in  the  cost  of  fossil  fuels,  alternate 
energy  sources  may  help  to  increase  the  profitability  of  greenhouses 
and  the  income  of  farmers  [4].  The  basic  strategy  of  the  greenhouse 
passive  heating  system  is  to  reduce  the  heat  losses  while  transferring 
the  excess  heat  from  inside  the  greenhouse  during  the  day  to  heat 
storage.  This  heat  is  used  during  the  night  to  satisfy  the  heating 
needs  of  the  greenhouse  [5,6]. 

Alternate  sources  of  energy  like  solar  energy  have  been 
proposed  for  heating  greenhouses  and  providing  the  optimum 
inside  conditions  during  winter  months.  Solar  energy  is  a  uni¬ 
versally  available  source,  but  its  practical  use  still  presents  some 
technical,  and  most  importantly,  economic  problems  [7  .  Some 
solar  systems  have  been  developed  but  they  are  too  expensive  for 
commercial  use  [7].  Several  types  of  passive  solar  systems  and 
techniques  have  been  proposed  and  used  by  [8—11  .  Bargach  et  al. 
[10]  used  a  solar  fat-plate  collector  to  improve  the  internal 
greenhouse  microclimate  and  tested  it  under  two  types  of  days; 
a  clear  day  and  a  cloudy  one.  Du  et  al.  12]  developed  a  model  to 
predict  the  performance  of  the  greenhouse  that  is  heated  with  a 
heat-pipe  system.  The  model  is  validated  with  experimental  data. 

Thermal  energy  can  be  stored  as  sensible  heat,  latent  heat,  reaction 
heat  or  combination  of  those  forms.  Many  thermal  energy  storage 
systems  and  applications  are  cited  in  literature  [13-17  .  The  most 
important  existing  greenhouse  heating  systems  that  use  a  sensible 
stored  energy  are:  water  storage  [18],  rock  bed  storage  [19,20]  and 
ground  storage  [21,22  .  The  latent  thermal  storage  has  many  advan¬ 
tages  over  the  sensible  one  such  as:  high  heat  capacity,  less  volume, 
low  storage  temperature,  thermal  energy  stored  and  released  at  an 
almost  constant  temperature.  There  is  a  large  number  of  Phase  Change 
Materials  (PCMs)  that  melt  and  solidify  at  a  wide  range  of  tempera¬ 
tures  which  are  used  in  greenhouse,  such  as  CaCl2  •  6H20, 
Na2SO410H2O  and  paraffin  [23-27].  Kern  and  Aldrich  [28]  utilized 
a  CaCl2  •  6H20  in  aerosol  cans  that  was  used  to  investigate  energy 
storage  possibilities  both  inside  and  outside  a  36  m2  ground  area  of 
the  greenhouse.  The  energy  storage  unit  inside  the  greenhouse 
(contained  1650  kg  PCM)  absorbed  the  excess  energy  of  warm  air 
inside  the  greenhouse  during  the  daytime.  The  energy  storage  unit 
outside  the  greenhouse  contained  1376.4  kg  PCM  and  two  solar  air 
collectors  with  8.55  m2  surface  area  each.  Results  of  this  study  showed 
that,  the  energy  stored  by  the  outside  unit  was  between  105.5  and 
158.25  MJ,  these  were,  21.1  and  31.65  MJ  for  the  internal  unit.  Huang 
and  Toksoy  [29]  placed  3000  kg  of  CaCl2  •  6H20  in  the  heat  exchanger 
inside  a  200  m2  glass  covered  greenhouse,  and  found  that  the  interior 
air  temperature  was  maintained  at  2  °C  higher  than  the  outdoor 
conditions.  Levav  and  Zamir  30]  used  CaCl2  •  6H20  in  greenhouses  and 


reported  that  the  required  air  temperature  in  the  greenhouse  was 
achieved  without  any  increase  in  the  relative  humidity.  Boulard  and 
Bailie  [31  ]  utilized,  in  a  greenhouse  with  176  m2  ground  areas,  double 
polycarbonate-cover  and  forced  ventilation,  2970  kg  of  CaCl26H20. 
The  PCM  was  packed  in  containers  and  placed  along  the  north  wall. 
The  PCM  could  provide  30%  night  heating  needs  during  the  winter 
period.  Boulard  et  al.  [32]  used  a  CaCl2  •  6H20  latent  heat  storage 
system  in  a  double-skin  polycarbonate  greenhouse  equipped  with  a 
forced  ventilation  system.  They  found  the  inside  greenhouse  tempera¬ 
ture  roughly  10  °C  higher  than  the  outside  temperature  during  typical 
nights  of  March  and  April.  Ozturk  et  al.  [33]  performed  an  experi¬ 
mental  evaluation  of  energy  and  exergy  efficiency  of  a  seasonal  latent 
heat  storage  system  for  greenhouse  heating  using  6000  kg  of  paraffin 
wax  as  a  PCM  in  a  180  m2  greenhouse  ground  areas.  They  observed 
that  the  average  energy  and  exergy  efficiency  of  the  arrangement  was 
40.4%  and  4.2%,  respectively. 

The  aim  of  this  work  is  to  determine  the  provided  heat  by  a 
Solar  Air  Heater  with  Latent  Storage  Collector  (SAHLSC)  inside  an 
Insulated  Greenhouse  with  Latent  Heat  System  (IGLHS).  Experi¬ 
ments  were  carried  out  in  two  greenhouses  designed  and  realized 
in  the  Research  and  Technology  Center  of  Energy  (CRTEn)  in 
Tunisia.  We  will  present  in  Section  2,  the  potential  of  solar  energy 
in  Tunisia.  In  Section  3,  we  will  describe  the  design  of  the 
experimental  greenhouses,  the  experimental  setup  and  measuring 
equipments.  In  Sections  4  and  5,  a  thermal  analysis  of  the  climate 
greenhouse  and  the  energy  balance  of  the  SAHLSC  are  carried  out. 
The  thermal  load  leveling  (TLL)  are  defined  in  Section  6.  The 
experimental  results  are  reported  in  the  Section  7.  Finally,  we  will 
present  the  economic  analysis  of  the  SAHLSC  in  Section  8.  The 
main  remarks  of  this  work  will  be  reported  in  the  conclusion. 


2.  Potential  of  solar  energy  in  Tunisia 

The  average  amount  of  solar  energy  radiating  in  Africa  is  about 
6  l<Wh/m2/day,  depending  on  the  latitude  of  the  country  [34].  In 
the  map  of  average  global  horizontal  solar  radiation  (Fig.  1),  we 
can  see  that  there’s  an  important  solar  potential  in  North-Africa. 

Tunisia  is  one  of  the  sunlight  countries,  it  lies  on  the  Medi¬ 
terranean  coast  of  Africa;  it  is  bordered  by  Algeria  to  the  west  and 
Libya  to  the  southeast;  with  the  following  coordinates:  Latitude 
36°43;N  and  Longitude  lO^'E.  Tunisia's  climate  is  temperate  in 
the  north,  with  mild  rainy  winters  and  hot  summers.  Tempera¬ 
tures  in  July  and  August  can  exceed  40  °C,  while  winters  are  mild 
with  temperatures  rarely  exceeding  above  20  °C  [35  . 

The  meteorological  station  on  the  site  of  Borj  Cedria  in  Tunisia 
permits  the  measurement  of  the  global  sun  flux  on  a  horizontal  plan, 
ambient  temperature,  wind  speed  and  wind  direction.  Figs.  2  and  3 
give  the  variation  of  meteorological  data  of  the  average  monthly  global 
solar  radiation  and  ambient  temperatures  in  a  period  of  2009-2011. 
We  can  see  that  the  global  solar  flux  has  the  same  trend  of  variation  as 
the  ambient  temperature.  It  is  observed  that  the  peak  of  average 
monthly  radiation  is  registered  in  July  (225  kW/m2)  with  a  highest 
monthly  ambient  temperature  surrounding  30  °C.  The  least  average 
monthly  radiation  is  65  kW/m2  with  a  monthly  ambient  temperature 
around  13  °C  in  December. 
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Nomenclature 

<7 

Stefan-Boltzmann  constant,  5.670  x  10  8  W/m2  K4 

T 

transmissivity 

A 

surface  area  (m2) 

Cp 

specific  heat  of  air  at  constant  pressure  (J/ kgK) 

Subscripts 

da 

density  of  air  (kg/m3) 

H 

heat  loss  coefficient  of  SAHLSC  (W/m2oC) 

A 

absorbed 

h 

heat  transfer  coefficient  (W/m2  °C) 

a 

ambient 

It 

total  solar  radiation  (W/m2) 

av 

average 

L 

latent  heat  (J/ kg) 

c 

cover  of  greenhouse 

LAI 

leaf  area  index 

ca 

canopy 

Lgh 

characteristic  length  of  greenhouse  (m) 

b 

bottom 

u 

collector  length  (m) 

ch 

charging 

l2 

collector  width  (m) 

dis 

discharging 

L3 

collector  depth  (m) 

e 

edge 

l 

characteristic  length  of  the  leaf  canopy  (m) 

finjeh 

final  of  the  charging  process 

lea 

average  length  of  the  canopy  (m) 

fin_dis 

final  of  the  discharging  process 

m 

mass  (kg) 

g 

glasses  of  SAHLSC 

m 

mass  flow  rate  (kg/s) 

H 

heat  gain 

Nu 

Nusselt  number 

i 

inside  greenhouse 

P(T ) 

saturated  water  vapor  pressure  at  temperature  T  (kPa) 

in 

inlet  of  SAHLSC 

Pr 

Prandtl  number 

inf 

infiltration 

Q 

heat  rate  (W) 

L 

liquid  phase 

Ra 

Rayleigh  number 

los 

leakage,  losses 

RE 

rate  of  air  infiltration  (m3/s) 

m 

melting 

Re 

Reynolds  number 

max 

maximum 

ra 

aerodynamic  resistance  (s/m) 

min 

minimum 

rs 

stomatal  resistance  (s/m) 

0 

outsider  greenhouse 

T 

temperature  (K) 

out 

outlet  of  SAHLSC 

t 

time  (s) 

PCM 

phase  change  material 

U 

velocity  (m/s) 

P 

packed-bed  absorber 

V 

volume  (m3) 

S 

solid  phase 

z 

depth  of  soil  (m) 

s 

soil 

sky 

sky 

Greek  symbols 

St 

storage 

t 

top 

a 

absorbtivity  for  solar  radiations 

u 

useful 

at 

absorbtivity  for  thermal  radiations 

V 

vertical  section  of  the  greenhouse 

Sb 

bottom  insulation  thickness  (m) 

Se 

edge  insulation  thickness  (m) 

Exposants 

A 

difference  in  time 

£ 

emissivity 

A 

absorbed  heat 

7 

psychometric  constant,  0.0667  kPa /K 

C 

convective  heat 

A 

thermal  conductivity  (W/m  K) 

Cd 

conductive  heat 

M 

thermal  conductivity  of  insulation  (W/m  K) 

L 

latent  heat 

n 

PCM  dynamic  viscosity  (Ns/m2) 

R 

radiation  heat 

Ps 

PCM  dynamic  viscosity  at  T  =  TPCm  (Ns/m2) 

p 

reflectivity 

3.  Materials  and  methods 

3.1.  Experimental  greenhouses  description 

In  this  study,  experiments  were  carried  out  in  two  small 
chapel-shaped  greenhouses  with  same  dimensions.  The  experi¬ 
mental  greenhouses  occupy  a  floor  area  equal  to  14.8  m2,  3.7  m 
wide,  4  m  long  and  3  m  high  at  the  center.  The  structure  is  all 
galvanized  steel  fixed  to  the  ground  with  stones  and  concrete. 
The  greenhouses  are  equipped  with  a  centrifugal  fan  controlled  by 
a  differential  thermostat.  The  fan  operated  when  the  temperature 
inside  the  greenhouse  exceeds  the  optimal  growth  temperature  of 
plant.  A  tomato  crop  (local  variety)  was  planted  in  both  green¬ 
houses  and  planted  at  7  plants  in  the  row  and  arranged  in  4  rows 
with  40  cm  between  the  rows  and  35  cm  between  the  plants. 


Physical  characteristics  of  the  greenhouse  glass  cover,  canopy  and 
soil  are  listed  in  Table  1.  The  reflectivity  and  transmissivity  spectra  of 
the  canopy  and  cover  layer  recorded  in  the  wavelength  range  (380- 
MOO  nm)  were  plotted  in  Figs.  4  and  5.  respectively.  The  reflectivity 
and  transmissivity  spectra  were  realized  using  LAMBDA950  UV-vis- 
NIR  spectrometer  equipped  with  an  integrating  sphere. 

The  first  greenhouse  is  named  a  Conventional  Greenhouse  (CG) 
where  its  walls  and  roofs  are  covered  with  a  3  mm  of  thick  glass 
(Fig.  6).  The  second  used  a  heating  system  and  was  named  as  an 
Insulated  Greenhouse  with  Latent  Heat  System  (IGLHS)  (Fig.  7(a)). 
The  south  oriented  wall  and  roof  of  the  IGLHS  are  covered  with  a 
3  mm  of  thick  glass.  Sidewalls  and  the  northern  roof  are  built  with 
0.4  m  and  0.6  m  of  thick  sandwich  panels  respectively. 

The  heating  system  of  the  IGLHS  is  a  new  Solar  Air  Heater  with 
Latent  Heat  Storage  Collector  (SAHLSC)  using  a  packed  bed  of 
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spherical  capsules  as  a  latent  heat  storage  system  [36].  A  new  type 
of  the  SAHLSC  was  designed  as  a  means  to  heat  the  interior 
environment  of  the  greenhouses  during  the  nighttime.  Fig.  7 
(b)  shows  the  SAHLSC  manufactured  inside  the  IGLHS.  The  length, 
the  width  and  the  total  volume  of  the  collector  are  2  m,  1  m  and 
0.28  m3,  respectively.  A  0.004  m  of  thick  transparent  glass 
cover  was  placed  0.015  m  apart  from  the  packed  bed  absorber. 
A  0.05  m  thick  polyurethane  insulation,  with  heat  conductivity  of 
0.028  Wm"1 1<~  \  is  placed  in  the  bottom  and  the  edge  of  the 
collector.  The  packed  bed  absorber  is  the  most  important  compo¬ 
nent  of  the  solar  air  heater  collector,  formed  of  spherical  capsules 
with  a  black  coating  and  fixed  with  steel  matrix.  Capsules, 
confined  with  a  PCM,  have  an  outer  diameter  of  0.077  m  and  are 
blow  molded  from  a  blend  of  polyolefin  with  an  average  thickness 
of  0.002  m  [36]. 

3.2.  Experimental  setup  and  measuring  equipments 

Experimental  greenhouses  (east-west  orientation),  are  con¬ 
structed  and  located  in  the  Research  and  Technology  Center  of 
Energy  in  Borj  Cedria,  in  Tunisia. 

During  the  sunshine  time  (charging  process)  a  fraction  of  the  total 
solar  radiation  received  inside  the  greenhouse  is  absorbed  by  the 
black  packed  bed  absorber  surface  of  the  SAHLSC.  The  absorbed 
thermal  energy  is  stored  as  a  sensible  and  a  latent  heat  forms  into 


the  collector.  At  night,  a  fan  blows  air  across  the  PCM  capsules 
and  extracted  the  stored  heat  into  the  IGLHS,  the  discharging 
process  done. 

K-type  thermocouples  inserted  inside  and  outside  the  greenhouses 
for  measuring  the  horizontal  distributions  of  temperature,  the  vertical 


Month 

Fig.  2.  Monthly  variation  of  the  global  solar  radiation  in  Tunisia  (2009-2011). 
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Fig.  1.  Global  solar  radiation  map  of  Africa  [34]. 
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Fig.  3.  Monthly  variation  of  the  ambient  temperature  in  Tunisia  (2009-2011). 


Table  1 

Cover,  canopy  and  soil  characteristic. 


Cover 

Canopy 

Soil 

Reflectivity  for  solar  radiations 

0.10 

0.1 

0.2 

Reflectivity  for  thermal  radiations 

0.10 

0.55 

0.16 

Emissivity 

0.9 

0.95 

0.9 

Transmissivity  for  solar  radiations 

0.85 

0 

0 

Transmissivity  for  thermal  radiations 

0.88 

0.4 

0 

Wavelength  (nm) 

Fig.  4.  Reflectivity  spectra  of  canopy  and  cover. 

distributions  and  the  temperatures  at  the  top  and  the  bottom  of  glass 
and  sandwich  panel.  The  temperature  and  relative  humidity  inside 
greenhouses  (HMP155A,  Campbell  Scientific  Inc)  were  measured  in 
the  center  at  a  height  of  1.5  m  from  the  ground.  The  HMP155A  sensor 
was  instaled  inside  a  14-plate  radiation  shield  houses  and  protects.  Its 
louvered  construction  allows  air  to  pass  freely  through  the  shield, 
thereby  keeping  the  probe  near  ambient  temperature  [37  .  Tempera¬ 
tures  under  the  soil  inside  the  greenhouse  at  0,  0.25  and  0.5  m  depth 
were  measured  using  three  PT-107  sensors.  The  average  temperature 
of  the  canopy  was  measured  using  a  non-contact  means  of  measuring 
the  surface  temperature  by  an  infrared  temperature  sensor  IR120.  The 
IR-SS  Solar  Shield  is  used  to  protect  the  IR120  sensor  from  direct  solar 
radiation  [37].  The  global  solar  irradiation  in  the  horizontal  plane  was 
measured  with  a  Kipp  and  Zonen  pyranometer  and  situated  at  1.5  m 


Wavelength  (nm) 


Fig.  5.  Transmissivity  spectra  of  canopy  and  cover. 


Fig.  6.  External  view  of  conventional  greenhouse. 


above  the  ground  in  the  center  of  the  IGLHS.  Ambient  temperature 
and  relative  humidity  are  measured  by  a  HMP155A  sensor  situated 
at  a  height  of  1.5  m  above  the  greenhouse.  Temperatures  along  the 
SAHLSC  were  measured  using  K-type  thermocouples.  Four  sensors 
were  inserted  inside  the  capsules  along  with  the  absorber  packed  bed. 
Also  four  thermocouples  are  used  to  measure  the  absorber  surface  of 
one  capsule,  the  glass,  the  inlet  and  the  outlet  air  temperatures  of  the 
SAHLSC.  All  climatic  parameters  and  measuring  instruments  were 
sample  recorded  every  10  min  using  a  CR5000  data  logger  (Campbell 
Scientific  Inc). 

The  external  global  solar  radiation,  wind  speed  and  direction 
are  provided  by  the  NRG  weather  station  in  the  Research  and 
Technology  Center  of  Energy  in  Tunisia.  This  station  is  equipped 
with  an  acquisition  system  that  record,  every  10  min  the  average, 
the  maximum,  the  minimum  and  the  standard  deviation  values  for 
each  sensor.  The  weather  station  provided  a  wind  speed  measure¬ 
ment  at  20  m.  The  wind  velocity  at  the  altitude  of  2  m  above  the 
greenhouse  was  corrected  by  using  the  power  law  equation  [38] 

(U,/U2)  =  (h,/h2y  (1) 

where  Hi  and  U2  are  the  wind  speeds  at  heights  /ii  and  h2 
respectively.  (3  is  the  power  law  exponent  and  the  subscripts  (1 
and  2)  correspond  to  different  heights. 

Uncertainty  analysis  is  needed  to  prove  the  accuracy  of  the 
experiments.  The  independent  parameters  measured  in  the 
experiments  reported  here  to  temperature,  air  velocity,  and  solar 
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Fig.  7.  (a)  External  view  of  IGLHS  and  (b)  SAHLSC  inside  greenhouse. 


radiation.  To  carry  out  these  experiments,  the  sensitiveness  of  data 
acquisition  system  is  about  +  0.001  °C,  the  measurement  error 
is  +  0.002  °C,  the  sensitiveness  of  the  thermocouple  is  +  0.01  °C, 
the  sensitiveness  of  the  IR120  is  ±  0.0004  °C  and  sensitiveness  of 
the  PT107  is  ±  0.01  °C.  The  HMP155A  errors  are  ±  0.02  °C  about 
temperature  and  +  2%  about  humidity.  An  anemometer  with 
+  0.01  ms-1  accuracy,  and  Kipp  and  Zonen  pyranometer  with 
+  3%  measurements  uncertainties  were  used.  The  sensitiveness 
was  obtained  from  a  catalog  of  the  instruments. 

4.  Thermal  analysis  of  greenhouse 

To  describe  the  thermal  behavior  of  the  greenhouse,  there  are 
four  components  that  play  important  roles  in  the  thermal  balance: 
the  cover,  the  canopy,  the  soil  and  the  inside  air  are  analyzed.  The 
following  assumptions  are  made: 

-  Radiation  heat  exchange  between  the  greenhouse  walls  and 
roofs  is  neglected. 

-  Storage  capacity  of  the  sandwich  panel  walls  and  roof  material 
is  neglected. 

-  Absorbtivity  and  heat  capacity  of  the  enclosed  air  is  neglected. 

-  Conduction  heat  exchange  between  air  and  the  sandwich  panel 
walls  and  roof  of  the  greenhouse  is  neglected. 

4.1.  Energy  balance  of  greenhouse  cover 

The  energy  balance  of  the  cover  of  greenhouse  described  by  the 
following  equation: 

Qc  +Qca-c  +  ^s-c“Q-c-o-Q.c-2  +  Q-c-s/<3/+Q.s-c  +  Q-ca-c  =  0  (2) 

where  Q hc  is  the  solar  radiations  absorbed  by  the  glasses  of  the 
greenhouse;  Q^a_c  is  the  solar  radiations  reflected  by  the  canopy 
and  absorbed  by  cover;  Q^-c  is  the  reflected  solar  radiations  by 
the  soil  surface  and  absorbed  by  cover;  Q^_0  is  the  convective 
heat  transfer  from  external  cover  of  the  greenhouse  to  ambient 
air;  o£_f  is  the  convective  heat  transfer  from  internal  cover  of  the 
greenhouse  to  inside  air;  Q^_sky  is  the  thermal  radiation 
exchanged  between  cover  and  sky;  Q*_c  is  the  thermal  radiation 
emitted  by  soil  surface  and  absorbed  by  cover  and  Q*a_c  is  the 
thermal  radiation  emitted  by  canopy  and  absorbed  by  cover. 

The  energy  balance  of  the  cover  Eq.  (2)  can  be  written  as 

ITAcac  +  IncUcPcaAca  +  h TCacpsAs  - hCc_ 0AC(T c  - T0) - hCc_ {AC{ Tc-T )) 

+  OecAc(TAsky  -Tc)-\-As(Jat,c£sTs4+Aca(Jat,c£caTca  =  0  (3) 

Different  relations  have  been  used  to  calculate  the  external 
heat  loss  coefficient  at  the  outside  surface  of  greenhouse  covers 
[39-42].  It  depends  on  the  greenhouse  geometry  under  any 
environmental  conditions.  Mac  Adams  [40]  suggested  a  relation, 


which  is  used  in  this  study 

hcc_0  =  5.7+3.8U0  (4) 

where  U0  is  the  wind  speed  outside  the  greenhouse. 

The  evaluation  of  internal  convective  heat  transfer  coefficient 
follows  Seginer  and  Livne  [42]  which  considers  the  combined 
effect  of  free  and  forced  convection. 

/l+  =  1.52|Tc-T0|1'/3  +  5.2  (+£)  ’  ^  (5) 

Uav  =  RE/Av  is  the  average  speed  of  the  air  inside  the  green¬ 
house;  RE  and  Av  are  the  rate  of  air  infiltration  and  surface  area  of 
the  vertical  section  of  the  greenhouse  respectively. 

4.2.  Energy  balance  of  canopy  surface 

The  energy  balance  of  canopy  surface  inside  the  greenhouse  is 
described  by  the  following  equation: 

0A  1  nR  _nR  —Qc  -Ol  =0  (6) 

^ca  '  -  ca  ca  *<- ca-i  ca  W 

where  Q^a  is  the  shortwave  radiation  absorbed  by  canopy;  Qj?_ca  is 
the  thermal  radiations  emitted  by  cover  and  absorbed  by  canopy;  Q*a 
is  the  thermal  radiation  emitted  by  canopy;  Q^a-i  is  the  convective 
heat  transfer  from  canopy  to  inside  air  and  Q^a  is  the  latent  heat 
transfer  between  the  canopy  and  the  inside  air  of  greenhouse. 

The  energy  balance  for  canopy  in  the  greenhouse  Eq.  (6)  can  be 
written  as 

T  ch^caAca  +  cco  T cAcat  ,ca  ~  £ca  &  T caACa  ~  bca  _  i Aca  ( T  ca  ~  T  i  ) 
daCpIAI(P\Tca)-P(Ti))A^  =  0 
Y  ra  +  rs  ca 

The  evaluation  of  the  convective  heat  transfer  coefficient  between 
the  inside  air  and  the  canopy,  Stanghellini  [43]  proposed 

hca  -  i  =  ~ ~ 1  /2  (Ca  I Tca  -  T,  |  +  207 1&)1  /4  (8) 

where  Uca  is  the  air  velocity  through  the  canopy  inside  the 
greenhouse. 

Uca  =  Uav(++)2/3  (9) 

The  saturated  water  vapor  pressure  P*(Tca )  at  the  temperature 
of  the  canopy  is  calculated  by  Tetens  equation  (Eq.  (10)),  the 
following  empirical  formula  relating  accurately  saturated  water 
vapor  pressure  to  temperature  for  the  temperature  range  between 

0  and  60  °C  [44,45] 

P*(Tt.)»0.6108exp(5^g|3)  „0) 
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The  aerodynamic  resistance,  ra  (sm-1)  of  the  canopy,  mainly 
depends  on  the  aerodynamic  regime  that  prevails  in  the  green¬ 
house  and  the  leaf  length.  It  is  expressed  as  [45] 


rn  —  220 


0.2 


U'i 


0.8 


01) 


where  l  is  the  characteristic  length  of  the  leaf  and  If,  is  the  mean 
interior  air  speeds. 

The  stomatal  resistance  rs  (sm_1 )  of  the  canopy  is  derived  from 
a  simple  empirical  relationship,  Eq.  (12),  with  global  radiation 
(thermal  and  humidity  dependences  of  greenhouse  plants  tran¬ 
spiration  were  neglected)  [45,46] 

rs  =  200  ^1  +exp(0  05(Tc/r  _ 50))  02) 


where  QA,  Qu,  Qst,  and  Q/os  are  the  absorbed,  useful,  stored  and 
lost  energies,  respectively. 

Based  on  Duffie  and  Beckman  [47],  the  useful  heat  gain  from 
the  collector  is 

Qu  =  mCp(Tout-Tin)  (19) 

The  radiation  absorbed  flux  of  the  absorber  packed-bed  is 
defined  as 

Qa  =  Ap(agTg)lT  (20) 

The  stored  heat  flux  during  the  charging  and  discharging 
phases  is  given  by 

Qch  =  [mpCMCpsiT  m  —T'ini_Ch.PCM)  +  mPCM^ 

+  mPCMCp,L(Tfin_ch,PCM  —  Tm)]/ At ch  (21) 


4.3.  Energy  balance  of  soil  surface 

The  energy  balance  of  soil  surface  inside  the  greenhouse  is 
described  by  the  following  equation: 

Q+Q+-Q"-Q+Q+  =  o  03) 

where  Q £  is  the  absorbed  solar  radiation  by  soil;  is  the 

convective  heat  transfer  from  soil  surface  to  air  inside  the  green¬ 
house;  Q^d  is  the  conductive  heat  transfer  into  soil  surface;  Q [R  is 
the  thermal  heat  radiation  emitted  by  soil  and  QR_S  is  the  thermal 
heat  radiation  emitted  by  greenhouse  cover  and  absorbed  by  soil. 

The  energy  balance  of  soil  in  the  greenhouse  (Eq.  (13))  can  be 
written  as 

1ttcolsAs  —hs_  i AS(TS  — T[ )  —  —AS(TS  —  Tz)  —  £scfTsAs  +  £coP  cAcols  =  0 

(14) 

The  convective  heat  transfer  coefficients  between  the  air  inside 
the  greenhouse  and  soil  surface  given  by  Seginer  and  Livres  [42] 

+  j  =  1 ,52|Ti  -  Ts  |  +  5.2  (j22J  (15) 


Qdis  =  [mPCM^p,L(Tm  —T fin_dis,PCM)  +  mPCM^ 

+  mPCM^p,s(Tin_dis,PCM  —  Tm)]/  A  tdis  (22) 

The  lost  heat  flux  is  given  by  Eq.  (23),  Hios  is  the  collector 
overall  heat  loss  coefficient.  The  thermal  energy  is  lost  from  the 
collector  to  the  surroundings  by  conduction,  convection  and 
infrared  radiation.  H/os  is  equal  to  the  sum  of  energy  loss  through 
the  top  (Hf),  bottom  (Hb)  and  edges  (He)  of  the  collectors  given 
below  [47] 

Qj05  =  HtosAc(h-Ti)  (23) 

Hios  =  Ht  3-Hb  +He  (24) 

The  top  loss  coefficient  from  the  collector  to  the  ambient  is 

=  (l/(/ip_g  +  /ip_g)  +  l/hg_(1)_1  (25) 

The  flow  is  laminar  during  the  charging  phase;  and  the 
appropriate  correlation  is  given  by  Churchill  [48] 

Nuch  =  2  +  (0.589  Pa1/4)/(l  +((0.469)/Pr)9/16)4/9 
Pr  >  0.7 


4.4.  Energy  balance  of  air  inside  greenhouse 

The  energy  balance  of  inside  air  of  greenhouse,  only  the 
convective  heat  transfer  is  considered.  It  is  assumed  that  no 
thermal  or  solar  radiation  is  absorbed  by  the  air.  The  energy 
balance  is  described  by  the  following  equation: 

Q.c-i  +  Qca-i  +  ^-1-  Qinf  +  Qh  =  0  (16) 

where  Q.^_,  is  the  convective  heat  transfer  from  inside  air  to  cover 
of  the  greenhouse,  Q£a_f  is  the  convective  heat  transfer  from 
inside  air  to  canopy  surface,  is  the  convective  heat  transfer 
from  inside  air  to  soil  surface,  Qfnf  is  the  heat  loss  due  to 
ventilation  or  infiltration  and  QH  is  the  heat  input  into  greenhouse. 

The  sensible  heat  balance  of  the  inside  air  of  greenhouse 
(Eq.  (16))  can  be  written  as 

hc_iAc(Tc  —  Ti)  +  hca  _  iACa(TCa  —  T{)  +  hs  _  iAs(T s  —  T,) 

-daCpRE(Ti-T0)  +  (lH  =  0  (17) 

where  QH  is  the  heat  gain  from  the  solar  air  heater  with  latent 
storage  collector  into  the  IGLHS. 


Then  the  flow  is  turbulent  during  the  discharging  phase;  the 
appropriate  correlation  is  given  by  Whitaker  44] 


Nudis  =  2  +  (0.4  fie1/2 +0.06  fie2/3)  Pr0A(g/gsf/4 
"0.71  <Pr<  380 
3.5  <  fie  <7.6  x  104 
1  <  <  3.2 


(27) 


hp_g  =  Nu(A/t) 

(28) 

hRP_g  =  ( a(T2p  +  Tg)(TP  +  Tg))/((1  /eP)  +  (1  /eg)  -  1 ) 

(29) 

hg -  Cl  =  £g(7(Tg  +  T^ky)(Tg  +  Tspy) 

(30) 

In  addition  the  energy  loss  through  the  bottom  (Hb)  and  the 
edges  (He)  is  as  follows: 

Hb  —  2j/ 5b 

(31) 

He  =  (Li  +L2)L3V(LiL2<5e) 

(32) 

5.  Energy  balance  of  the  SAHLSC 


6.  Thermal  load  leveling  (TLL) 


The  energy  analysis  balance  during  both  the  charging  and 
discharging  phases  of  the  SAHLSC  is 


(18) 


The  thermal  load  leveling  (TLL)  is  a  relative  index  to  represent 
the  fluctuation  of  temperature  inside  the  greenhouse.  The  less  the 
fluctuations,  the  better  is  the  environment  for  plants  inside  the 
greenhouse.  For  thermal  heating  of  the  greenhouse  with  minimum 


Qa  —  Qu  +  Qst  +  Qios 


38 


S.  Bouadila  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  35  (2014)  31-41 


fluctuations  the  TLL  should  has  a  minimum  value  [49]. 


TLL  = 


Ti,  max  Tj  m[n 
T i,  max  +  T i  m[n 


(33) 


7.  Results  and  discussions 

7.2.  Comparative  study  of  greenhouses  climate 

In  the  CG  and  IGLHS,  tomato  plants  are  transplanted  on  the 
23th  of  January  and  the  harvesting  started  on  the  15th  of  April 
2013.  The  major  environmental  factors  that  affect  the  growth  and 
the  precocity  of  the  production  of  greenhouse  plants  are  carbon 
dioxide,  light,  humidity  and  temperature.  The  tomato  plants  have 
two  optimum  temperatures,  one  during  the  day;  which  varied 
between  20  and  30  °C,  and  the  other  nocturnal  temperature  which 
is  the  most  crucial  temperature,  varied  between  12  and  18  °C. 

During  the  growing  period,  the  February  has  the  least  average 
monthly  radiation  is  85  kW/m2  with  a  monthly  ambient  temperature 
around  12.5  °C  (Figs.  2  and  3).  An  experimental  study  has  been  carried 
out  during  the  month  of  February  2013  which  is  the  winter  season. 

Fig.  8  presented  the  instantaneous  registered  air  temperature 
inside  CG  and  IGLHS  (lst-28th  February,  2013).  The  fluctuation  of 
air  temperature  inside  the  two  greenhouses  is  closely  attributed  to 
weather  conditions  such  as  solar  radiation  and  outside  tempera¬ 
ture.  Under  the  experimental  conditions,  the  nocturnal  variation  of 
temperatures  inside  greenhouses  varied  between  2  and  5  °C.  The 
method  of  utilizing  the  new  solar  air  heater  with  latent  storage 
collector  as  a  heating  system  inside  the  greenhouse  shows  mean¬ 
ingful  value  in  February. 

An  experimental  study  has  been  carried  out  on  typical  selected 
days  of  the  system  operation  from  23th  to  27th  February,  2013  in 
winter  season. 

The  external  global  solar  radiation  and  wind  speed  for  the  5- 
day  period  in  the  region  of  the  greenhouses  are  shown  in  Fig.  9. 
The  maximum  of  the  global  solar  irradiation  intensity  in  the 
horizontal  plane  varied  between  300  and  650  Wm~2.  The  two 
first  days  of  this  period  were  very  windy,  air  speed  varied  between 
3  and  8  ms-1.  The  average  speed  of  other  days  is  of  2  ms-1. 

The  daily  variation  of  the  ambient  temperature,  air  tempera¬ 
ture  inside  CG  and  air  temperature  inside  IGLHS  are  presented  in 
Fig.  10  (for  the  same  5-day  period).  The  air  temperature  inside  the 
greenhouses  followed  the  same  trend  as  the  ambient  temperature 
and  exceeds  the  daily  tomato  optimum  temperature  under  severe 
winter  weather.  Consequently,  it  is  necessary  to  improve  the  air 
conditioning  of  this  greenhouse.  A  new  solar  air  heater  collector 
with  a  latent  heat  storage  system  is  integrated  inside  the  IGLHS. 
This  collector  is  used  to  store  the  excess  of  the  transmitted  solar 
irradiation  in  the  packed  bed  absorber  and  to  provide  it  at  night. 
The  discharging  process  of  the  SAHLSC  starts  at  21  h  (local  time). 
A  fan  was  used  to  blow  the  air  at  a  fixed  speed  equal  to  lm/s. 

The  23th  and  27th  of  February,  2013  is  characterized  by  clear 
sky  with  maximum  ambient  temperature  equal  to  21  °C.  We 
observed  that  the  diurnal  temperature  attained  39  °C  and  41  °C 
respectively  into  CG  and  into  IGLHS.  The  nocturnal  variation  of 
temperature  inside  the  IGLHS  exceeds  the  temperature  inside  CG 
with  a  difference  of  5  °C  from  21  h  (local  time)  until  the  morning 
after  (Fig.  10).  On  27th  of  February,  2013  the  SAHLSC  maintained 
the  air  temperature  inside  IGLHS  constant,  around  12  °C,  along  the 
night,  due  to  the  fraction  of  solar  radiations  stored  in  the  packed 
bed  absorber.  However,  the  temperature  inside  CG  decreases 
progressively  to  the  ambient  temperature  after  2  h  from  sunset. 
The  24th  February,  2013  was  characterized  by  severe  solar  radia¬ 
tion  fluctuations,  the  ambient  air  temperature  does  not  exceed 
11  °C  and  wind  velocity  varied  from  4  to  8  m/s.  We  noted  that  the 
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Fig.  8.  Air  temperature  inside  CG  and  air  temperature  inside  IGLHS  as  a  function  of 
days.  (lst-28th  February,  2013). 
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Fig.  9.  Global  solar  radiation  and  wind  velocity  as  a  function  of  days  (23th-27th 
February,  2013). 
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Fig.  10.  Ambient  temperature,  air  temperature  inside  CG  and  air  temperature 
inside  IGLHS  as  a  function  of  days  (23th-27th  February,  2013). 


temperature  into  the  IGLHS  exceeds  the  temperature  inside  CG 
with  a  difference  of  3  °C  all  the  night. 

Fig.  11  presents  the  effect  of  the  SAHLSC  on  relative  humidity 
inside  the  greenhouse.  At  night,  the  IGLHS  permitted  to  maintain 
the  air  relative  humidity  to  10-20%  lower  than  the  CG.  The  SAHLSC 
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Fig.  11.  Humidity  inside  CG  and  IGLHS  as  a  function  of  days  (26th-27th 
February,  2013). 
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Fig.  13.  Nighttime  heating  recovered  from  SAHLSC  and  IGLHS  heating  require¬ 
ments  as  a  function  of  days  (26th-27th  February,  2013). 
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Fig.  12.  Useful,  absorbed  and  stored  heat  rates  as  a  function  of  day  times  inside  the 
IGLHS  (27th  February,  2013). 


creates  a  passive  dehumidification  process  at  night  time  due  to  the 
increase  of  the  air  temperature  inside  the  IGLHS. 

7.2.  Performance  of  the  SAHLSC  inside  IGLHS 

Variations  of  useful,  absorbed  and  stored  energy  are  presented 
in  Fig.  12  during  27th  of  February,  2013.  Therefore,  to  evaluate  the 
stored  energy,  the  energy  balance  equation  (Eq.  (18))  has  been 
used.  A  great  part  of  the  absorbed  solar  heat  was  stored  inside  the 
PCM.  It  is  observed  that  during  the  initial  period  of  charging,  the 
instantaneous  heat  stored  increases  with  insolation  and  towards  a 
maximum  value  of  0.7  kW  at  13:00,  we  noted  that  at  this  time  the 
absorbed  heat  is  1.3  kW,  53%  of  the  solar  energy  is  stored  into 
the  collector.  As  the  discharging  process  proceeds,  the  PCM 
starts  solidifying  and  the  used  heat  is  uniform  for  a  longer  period. 
The  uniform  value  of  the  used  heat  is  between  550  and  300  W  all 
the  night. 

Energy  balance  equation  (Eqs.  (2),  (6),  (13),  (17)  and  (18))  have 
been  used  for  calculating  requirements,  excess  and  recovered 
thermal  energy  (MJ)  inside  IGLHS.  Fig.  13  presented  the  daily 
heating  requirements  and  the  nighttime  heating  recovered  from 
SAHLSC,  during  the  5-days  period  (23th-27th  February,  2013).  The 
stored  energy  in  the  solar  air  heater  with  thermal  storage  is 
13.7  MJ  on  23th  February,  2013.  The  contribution  of  stored  energy 
in  the  SAHLSC  to  the  nighttime  heating  requirements  was  31%,  27% 
and  26%  respectively  on  23th,  26th  and  27th  February,  2013  under 
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Fig.  14.  Thermal  load  leveling  of  CG  and  IGLHS  as  a  function  of  days  (26th-27th 
February,  2013). 

best  ambient  conditions.  On  24th  and  25th  February,  2013,  the 
recovered  and  stored  heat  is  6%  and  11%  with  requirements  heat 
exceeds  65  MJ. 

The  performance  of  IGLHS  with  a  SAHLSC  heating  system  has 
been  evaluated  in  terms  of  thermal  load  leveling,  TLL,  using  Eq.  (3). 
The  TLL  is  used  to  quantify  the  fluctuations  of  temperature  inside 
the  greenhouse.  In  winter,  TLL  should  have  lower  value  by  incor¬ 
porating  a  heating  system  due  to  the  increase  of  (T,  max  +  ^i  min )  as 
well  as  decrease  of  (T,  max  -  T,  min )  as  compared  to  TLL  without 
heating  system.  The  results  of  the  daily  variation  of  thermal  load 
leveling  for  CG  and  IGLHS  greenhouse  have  been  shown  in  Fig.  14. 
It  can  be  seen  that  the  value  of  TLL  is  maximum  for  CG  greenhouse 
and  it  is  reduced  to  25%  for  IGLHS  on  27th  February,  2013.  The 
lower  values  of  thermal  load  leveling  indicate  the  decrease  in  the 
fluctuations  of  greenhouse  air  and  thereby,  there  occurs  an 
improvement  of  desired  environment  for  plants  in  the  greenhouse. 


8.  Economic  analysis 

The  economic  analysis  (SAHLSC  cost,  energy  cost  and  payback 
period  in  comparison  with  the  expected  service  life  cycle  of  the 
installation)  is  done  without  taking  into  account  the  interest  and 
the  inflation,  we  suppose  that  all  the  cost  is  paid  without  credit 
payment  [50  .  The  life  cycle  of  the  SAHLSC  is  considered  20  years 
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and  its  thermal  performance  degradation  is  assumed  to  be  1% 
per  year. 

In  order  to  evaluate  the  economic  analysis,  the  SAHLSC  system 
is  compared  to  non-solar  alternative  heating  system  (fuel  boiler). 
In  this  study,  two  systems  are  considered: 

-  Hybrid  system:  the  SAHLSC  and  the  fuel  boiler.  The  fuel  boiler 
is  the  auxiliary  system  meeting  the  energy  needs. 

-  Conventional  system:  the  most  heating  system  used  for  green¬ 
houses  conditioning  is  the  fuel  boiler. 


8.1.  The  total  cost  system 

The  total  cost  of  the  greenhouse  heating  system  includes  the 
investment  cost  and  the  energy  consumption  cost.  The  investment 
cost  includes  the  system  cost  and  the  operation  and  maintenance 
cost  [51,52]. 

The  total  cost  of  the  conventional  system  is  defined  by  adding 
the  initial  investment  cost  of  the  fuel  boiler  and  the  fuel  con¬ 
sumption.  For  the  hybrid  system  the  total  cost  is  the  sum  of  all  the 
initial  investment  cost  (SAHLSC + fuel  boiler)  and  the  fuel  con¬ 
sumption  without  the  supplied  heat  cost  by  the  SAHLSC. 

-  The  system  cost  is  defined  by  the  equipment,  the  realization 
and  the  installation  costs.  The  hybrid  system  cost  includes  the 
fuel  boiler  and  SAHLSC  costs.  Since  the  SAHLSC  is  not  yet 
commercialized,  there  is  no  official  price  for  this  system.  Thus, 
the  capital  cost  of  the  solar  collector  was  roughly  determined 
as  the  costs  spent  to  develop  the  experimental  system. 

-  The  operation  and  maintenance  cost  refers  to  heat  production 
coupled  with  full  load  hours.  The  electricity  consumption  for 
operating  the  fan  was  considered  as  operating  costs.  The 
maintenance  cost  of  the  SAHLSC  is  considered  negligible. 


8.2.  Energy  consumption  for  the  greenhouse  heating 

The  sizing  of  the  greenhouse  heating  system  is  based  on  the 
nocturnal  setpoint  temperature  (12  °C).  This  temperature  is  sui¬ 
table  for  tomato  cultivation  implemented  in  the  experimental 
greenhouse.  The  necessities  of  the  energy  to  maintain  the  noctur¬ 
nal  air  temperature  under  the  greenhouse  at  the  setpoint  value  is 
equal  to  the  energy  losses.  The  total  energy  losses  are  determined 
by  adding  the  energy  losses  of  the  greenhouse  components 
(cover,  ground,  infiltration,  and  canopy).  The  Tunisian  greenhouse 
heating  period  is  from  November  to  March,  equivalent  to  1500  h 
of  heating.  In  this  period,  the  fuel  boiler  in  the  hybrid  system, 
functions  only  600  h.  Therefore  it  ensures  an  economy  of 
1800  kWh/year. 

8.3.  The  payback  duration 

The  aim  of  the  economical  analysis  is  to  find  out  the  payback 
period  of  the  SAHLSC.  The  payback  duration  indicates  the  time 
that  the  SAHLSC  generates  the  amount  of  energy  equivalent  to  the 
energy  provided  during  the  life  cycle  [53  . 

The  cost  of  each  heating  systems  (hybrid  and  conventional 
system)  and  the  payback  time  of  the  SHALSC  is  mentioned  in 
Table  2.  As  a  result,  the  payback  period  of  solar  air  heater  with 
latent  storage  is  found  as  5  years. 

8.4.  Environmental  impact  of  the  SAHLSC 

From  2005  to  2008,  the  greenhouse  areas  of  Tunisia  increased 
from  4600  ha  to  9000  ha  [54],  which  consume  3  x  105  t  of  fuel 


Table  2 

Comparison  of  the  cost  and  the  payback  times  between  SAHLSC  system  and  the 
fuel  boiler. 


Hybrid  system 
(SAHLSC  +  fuel  boiler) 

Conventional  system 
(fuel  boiler) 

System  cost  ($) 

3600 

1800 

Operation  cost  ($) 

4.75 

116.55 

Maintenance  cost  ($) 

150 

150 

Annual  energy 

1200 

3000 

consumption  (kWh) 

Gain  of  energy  (kWh) 

1800 

- 

Energy  cost  ($) 

259.2 

651.25 

Total  cost 

4014 

2717.8 

Payback  of  SAHLSC 

5 

- 

(year) 

during  5  months  of  heating  period.  The  maximum  of  the  C02 
content  is  about  16%  of  the  fuel  consumption  [55],  therefore  the 
conventional  heating  system  causes  the  emission  of  3.12  x  105  t  of 
C02  which  is  equal  to  2.8  x  10-2  t  per  person  [56].  The  reduction 
of  carbon  dioxide  emissions  using  the  SAHLSC  may  be  a  promising 
solution.  Indeed,  the  SAHLSC  ensure  the  greenhouse  conditioning 
during  3  months,  it  can  reduce  the  C02  emission  at  1.68  x  10“ 2 1 
per  person  which  presents  1%  of  the  total  of  the  C02  emissions. 


9.  Conclusion 

An  experimental  study  was  conducted  to  evaluate  the  night¬ 
time  recovered  heat  of  the  SAHLSC  in  the  east-west  oriented 
greenhouse.  During  winter  months,  the  maximum  of  solar  radia¬ 
tion  falls  on  the  south  wall.  Therefore,  a  SAHLSC  is  proposed  in  this 
study  for  storing  the  excess  of  the  transmitted  solar  irradiation  in 
the  packed  bed  absorber  and  to  provide  it  at  night.  This  stored 
energy  is  extracted  to  the  greenhouse  by  convection  and  radiation 
heat.  The  performance  of  the  solar  system  was  evaluated  in  terms 
of  thermal  load  leveling  and  nighttime  recovered  heating  for 
a  5-day  period  of  February,  2013  in  Tunisia.  As  a  result  of  this 
system,  the  amount  of  the  nighttime  recovered  heat  of  this  system 
attains  31%  of  the  total  requirements  of  heating.  The  relative 
humidity  was  found  to  be  an  average  10-20%  lower  at  night  time 
inside  the  heated  greenhouse.  The  nocturnal  variation  of  tempera¬ 
ture  inside  the  IGLHS  exceeds  the  temperature  inside  the  CG  by 
5  °C.  The  solar  air  heater  remains  a  uniform  useful  heat  during  the 
nighttime.  The  uniform  value  of  the  used  heat  is  between  550  and 
300  W  all  the  night.  This  useful  heat  was  not  affected  by  the  severe 
global  solar  radiation  fluctuation  during  the  charging  mode.  The 
evaluation  of  the  economic  analysis  of  the  SAHLSC  assessed  that 
the  payback  period  was  approximately  5  years  if  the  system 
operates  only  3  months  per  year.  The  environmental  impact  from 
reducing  carbon  emissions  is  about  1%  of  the  total  Tunisian 
greenhouses  emissions. 
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